Superconductivity near Itinerant Ferromagnetic Quantum Criticality 
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Superconductivity mediated by spin fluctuations in weak and nearly ferromagnetic metals is 
studied close to the zero-temperature magnetic transition. We solve analytically the Eliashberg 
equations for p-wave pairing and obtain the normal state quasiparticle self-energy and the super- 
conducting transition temperature T c as a function of the distance to the quantum critical point. 
We show that the reduction of quasiparticle coherence and life-time due to scattering by quasistatic 
spin fluctuations is the dominant pair-breaking process, which leads to a rapid suppression of T c to 
a nonzero value near the quantum critical point. We point out the differences and the similarities of 
the problem to that of the theory of superconductivity in the presence of paramagnetic impurities. 
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Superconducting ground states have recently been dis- 
covered in materials sitting close to the zero-temperature 
phase boundary of both antiferromagnetic and ferromag- 
netic transitions |l],|j . It is quite natural to suspect that 
paring in these materials is mediated by spin-fluctuations 
which are enhanced near the magnetic quantum critical 
point (QCP) However, it is known that quasiparti- 

cles scattered by critical spin-fluctuations are subject to 
severe non-Fermi liquid self-energy corrections that are 
in general pair-breaking. The interplay between these 
two competing effects, generic to superconductivity near 
quantum phase transitions, has not been fully understood 
and is the subject of this work. 

We focus on the case of p-wave superconductivity near 
the ferromagnetic QCP where the Curie temperature is 
driven to zero by, e.g. pressure or doping [p|JlC[|. In the 
vicinity of the transition, the metals can be described by 
weak ferromagnetic (FM) and nearly FM Fermi liquids 
at low temperatures, respectively pljp^| . The impor- 
tant low energy excitations are the quasiparticles and the 
long wavelength spin fluctuations. Since the quasiparti- 
cles couple strongly to spin fluctuations near the QCP, 
it is necessary to adapt the framework of the strong cou- 
pling Eliashberg equations. The long wavelength nature 
of the critical spin fluctuations near the ferromagnetic 
QCP jl3| makes the solution of these equations theoret- 
ically tractable, enabling us to obtain analytical results 
for the self-energy and the transition temperature T c . We 
find the onset of p-wave superconductivity in both the 
Fermi liquid T c < T* and the quantum critical regime 
T c > T* , where T* is the characteristic frequency for 
spin fluctuations. We show that T c is rapidly reduced 
on approaching quantum criticality with T* — > due to 
the rapid reduction of quasiparticle coherence and life- 
time caused by quasistatic scattering of spin fluctuations 
analogous to the suppression of T c by paramagnetic im- 
purities H]. Interestingly, T c remains finite at the FM 
QCP. We shall not discuss the problem of coexistence 
of ferromagnetism and superconductivity here, because 



the latter requires a careful treatment of the coupling 
between the superconducting order parameter and the 
electromagnetic fields associated with the ferromagnetic 
fluctuations Q that is beyond the scope of this work. 

In a FM Fermi liquid with spin polarization along the 
z-axis, the single-particle Green's function has the form, 
G a {e,p) = a r7 /[e~v rT (p-p cr ) + if]sgn(e)}, where r\ -> 0+, 
p a (a =|, I) are the Fermi momenta of the spin up and 
down electrons, and a a are the wave-function renormal- 
izations. For a weak FM metal, such as the one close to a 
continuous FM transition, the difference in the Fermi mo- 
menta is small, S = |pj — pjj <C Pi, I' an< ^ ^ ^ s possible to 
set p-f = Pi — pf and a-f = Jll[]. The spin fluctuations 
are described by the propagators of the electron spin- 



density, Dij (x, t) 



-i{(Si%r)- {8^(0,0)- (Sj))), 



i,j = x, y, z. The long wavelength and low energy spin- 
fluctuations in the transverse and longitudinal channels 
are given by Q , 
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Here a w (S/pf) 2 <§C 1 measures the distance to the crit- 
ical point, S = vfSNf/2 is the averaged uniform spin 
density, A = 2vf Pf is an energy scale of the order of the 
Fermi energy, and Np — p 2 F / (2tt 2 vf) is the density of 
states per spin at the Fermi level. Notice that the trans- 
verse spin-wave only emerges for q -C 5 with a dispersion 
u) s (q) ~ v F d~(q/2p F ) 2 . 

On the paramagnetic (PM) side, the spin rotation 
symmetry is restored. The spin fluctuations become 
isotropic, D± = D\\ = D. All three modes take on the 
paramagnon form in Eqs. (Q) with a determined by the 
spin correlation length £ s , a ~ £~ . 

FM spin fluctuations are known to be pair-breaking 
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in the s-wave channel and therefore suppress the conven- 
tional phonon-mediated superconductivity in PM metals. 
This problem was studied by Berk and Schrieffer [jl6| us- 
ing the strong coupling Eliashberg theory and is believed 
to be the reason for why transition metals close to the 
FM instability have, if any, a low T c . In contrast, the 
pairing interaction is attractive in the / = odd angular 
momentum channel, raising the interesting possibility of 
FM spin fluctuation mediated spin-triplet superconduc- 
tivity JT^ ]. In general, the presence of a spontaneous 
magnetization in the FM phase leads to a set of four 
Eliashberg equations. However, for weak ferromagnets 
with small moments, the spin dependence of the self en- 
ergy and the gap function can be ignored. For notational 
convenience, we shall limit the presentation to the PM 
phase, the modification on the FM side is straightforward 
and the differences will be noted explicitly. 

Expressing the self energy in the Nambu formalism, 
E p (w) = [1 — Z p (u))]ustq + Z p ((jj)A p (uj)ti, where A p (uj) is 
the gap function, the linearized Eliashberg equations for 
the self-energy and the gap function are given by fig], 



xD(p — k, 
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Z p (ioj n )A p (iuj n ) = g 2 siT^ 



xD(p-k,iuj n -ie n ), (3) 
d 3 k Z k (ie n )A k (ie n ) 
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Here, go is the coupling constant of the quasiparticles to 
spin fluctuations and £ k = k 2 /2m—fi. Note that the pres- 
ence of Heisenberg symmetry in the PM phase guarantees 
three identical soft modes contributing to the self-energy 
but one longitudinal mode to the gap equation for triplet 
(p-wave) pairing ||, i.e. sq = 3 and s± = 1. To study 
the influence of the departure from this symmetry on T c , 
such as in the case of Ising spins with so = Si = 1, we 
keep So and si as general parameters On the FM side, 
spin rotation symmetry is broken. While the gap equa- 
tion (Q) only contains the longitudinal mode and stays 
invariant, the self-energy equation ([|) needs to be mod- 
ified in a straightforward manner to reflect the different 
contributions from the transverse modes. 

It is customary to proceed by separating the k- 
integral into angular and amplitude parts, J d 3 k — > 
{Ap 2 f /vf) J xdx J d(f> J d£ k , where x — sin(6>/2), 8 is the 
angle between k and q — p — k, and £ k satisfies the kine- 
matic constraint q 2 ~ Ap 2 F x 2 +^ k /v 2 F . As it will turn out 
later, close to the QCP and compared to their frequency 
dependence, the self energies have a negligibly weak mo- 
mentum dependence. We thus drop the /c-dependence in 
Z k (iuj) — > Z(iuj) and project the gap function into the l- 
th angular momentum channel, A k (iu) — > Ai(iui). With 
this approximation, it is possible to carry out the inte- 
gral over In order to make analytical progress, we 
analytically continue to real frequencies and obtain, 
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P l (l-2x 2 ), (6) 



where g 2 = g^Np/2, Pi(x) is the Legendre polynomial, 
and U 2 (e) = a + x 2 — ine/AxA. 

We first solve Eq. m) to derive the quasiparticle self 
energy in the normal state. The self consistency in the 
this equation is crucial near the QCP, since as a — > 0, 
the important low frequency cutoff in the denominators 
is the self energy itself. Carrying out the integrals, it is 
straightforward to show that the dominant contributions 
to the self energy come from the scattering by the spin 
fluctuations with momentum transfer q 3> S and having 
the same form in both the PM and the FM phase close 
to the QCP. We find that the characteristic energy scale 
for spin fluctuation, T* ~ a 3 / 2 A, enters as an important 
crossover temperature scale. For y = max(T, e) < T* , 
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the self energy behaves as in a Fermi liquid, 

S(e,T) « -c'eln(A/T*)-ic"y 2 /T*, y<T*, (7) 

where c',c" ~ sog 2 - However, for y > T*, the scatter- 
ing by spin fluctuations is enhanced and the self energy 
becomes non Fermi liquid like with the real part, 

£'(e,T)«-c'em(A/y), y > T*. (8) 

This leads to a quasiparticle residue Z, Z^ 1 = 1 — 
[<9£'/(9e]| e=0 = I + c'ln[A/Max(T*,T)], that vanishes log- 
arithmically on approaching the QCP as in a marginal 
Fermi liquid p9| . For e > T 3> T* , the imaginary part 
of £ follows, 
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E"(e) « -c"7re/2, e>T>T*. 



(9) 



In the quasistatic regime, T > e ^> T*, another energy 
scale arises by comparing the coherence length £ co h ~ 
vf/T to the spin correlation length £ s ~ 1/y/a. We find 
for £ co h > £,s, 



T T A T* 
E"(T)«-c"Tln — , £l n il«(i^)i/3. 



(10) 



For even smaller T*, such that £ co h < £ s , the quasipar- 
ticles scatter off essentially uncorrelated spins. In this 
case, the unphysical singularity in Eq. ([l0|) as T* — > 
must be removed by the self-consistency of Eq. (||) where 
E" itself becomes the cutoff. We find 



E" « -c"Tln(A 2 T/|E"| 3 ) + 3c"V^TA/2|E"| 
which has the following self-consistent solution, 



(11) 



E"«- C "Tln^ + |c"VSA, |laA >( ^)i/3. (i 2 ) 

The same behavior of E holds on the FM side since the 
dominant scattering process has q 5. In this range of 
q, the breaking of spin rotation invariance is insignificant 
as seen in Eqs. (Q) and (j^). The spin wave contribu- 
tion is not important since the difference between the 
Fermi momentum (S) is much larger than most of the 
momenta carried by the spin waves. Notice that the in- 
elastic scattering rate in quasistatic limit increases with 
a square-root singularity as a — > 0, leading eminently 
to a rapid suppression of the superconducting transition 
temperature T c on approaching the QCP. 

We next determine T c by solving Eqs. (^) for the sim- 
plest 1 = 1, p-wave case [^(J To treat the effects of 
both mass renormalization and scattering lifetime, we 
write Z = Z' + iZ" and the complex gap function as 
A = A' + iA" . Taking the imaginary part of the gap 
equation (^), we obtain to leading order in T c /To, To ~ A 
being the cutoff frequency for spin fluctuations - a mag- 
netic analog of Debye frequency, 



Z'(lu)A"(uj) + (3Z"(lu)A'(lu) Si /s ~ 0. 



(13) 



Here (3 — sq/si — 1 reflects the spin symmetry. Writ- 
ing for small u, uZ(u) = [1 + \(T)]ui — iT(oj), where 
A = — E'/o; is the effective coupling and T = — E" is 
half the inverse life-time of the quasiparticles, Eq. ( O ) 
becomes f3T(u))A'(u>) = (1 + X)ujA"sq/si, which allows 
us to account for the damping of the order parameter in 
terms of a real effective gap function, 



A eff (w) = A'(w)[l + ^/w) s 



(14) 



where T = (3siT/sq(1 + A). In contrast to A', A e gf has 
a weaker w-dependence and remains finite in the small 
uj limit. Now we can rewrite the real part of the gap 
equation ([j]) as an integral equation for A e g , 
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tanh(^) , eA ? ff(£) Re, ' 



'2T C ' e 2 + ri(e,T c ) U^{oj-t, X y 



(15) 



The scattering rate enters as a low-energy cutoff of the 
logarithmic singularity responsible for the superconduct- 
ing instability. Next we attempt an approximate analyti- 
cal solution of T c from Eq. (15) where the weak frequency 
dependence of A e ff can be neglected. 

Fermi liquid regime. Consider first the case T c <C 
T* -C Toj i-e. the onset of superconductivity in the Fermi 
liquid regime away from the QCP. Since the temperature 
is much lower than the characteristic spin fluctuation fre- 
quency, inelastic scattering dominates, but with the ordi- 
nary Fermi liquid scattering rate (see Eq. that is much 
smaller than max(/csT, e). The effect of a nonzero T on 
T c is thus small and negligible. Solving Eq. jl5| ) to next 
to leading order in T*/T , we obtain 



-[V+/3(2A+3)+A 2 ]/(2A-3) 



(16) 



where A = \nT /T*, rj' = (3/2s iff 2 ) + (tt 2 /24) + 
6\n(V2j/n) — 3, and In 7 ~ 0.577 is Euler's constant. 
As T* is reduced towards the QCP, spin fluctuations in- 
crease and pairing is enhanced. This causes T c to rise 
initially. Reducing T* further eventually causes the sys- 
tem near T c to lose sensitivity to the finite correlation 
length, leading to new physics associated with supercon- 
ductivity near quantum criticality. 

Quantum critical region. Here, Tq ^> T c ^> T* , the 
superconducting transition occurs inside the quantum 
critical regime. Since the temperature is much higher 
than the characteristic quantum spin fluctuation energy, 
inelastic scattering is negligible and the dominant pair- 
breaking effect comes from quasistatic (to < T c ) spin fluc- 
tuations with a scattering rate T(T) = —f3siY"(T)/s [l + 
A(T)]. The suppression of T c due to T is thus reminiscent 
of Abrikosov and Gor'kov's theory of superconducting 
alloys with paramagnetic impurities |14"| . Accordingly, 
Eq. ( |l5| ) has the solution, 



r(T e ) 

2vrT r 



(17) 



where ip is the digamma function and T c q is the transi- 
tion temperature in the absence of T. To leading order 
in Tc/Tq, 



^0 = -[/3 +v /^+2/3[l-^(^) 2 /3]+,+ ^(^. 



2/31 



Tn 



(18) 



with (3 = 3 + P and n = (3/2s ; .g 2 ) + 61n(27/7r) - tt 2 /24. 
At the QCP, we find T c0 /T ~ 10~ 5 for Heisenberg sym- 
metry and ~ 10 ~ 3 for Ising symmetry. 
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A few remarks are in order for T c in the quantum crit- 
ical regime, (i) For Ising spins, (3 = 0. Eq. @ shows 
that T c is not affected by a finite quasiparticle life-time, 
i.e. T c ~ T c q. Furthermore, to leading order, T c q is not 
reduced by the real part of the self energy. This is in fact 
a manifestation of Anderson's theorem for nonmag- 
netic impurities. It arises in our case from the cancella- 
tion of the self-energy effects in the gap equation ([l5]) in 
the quasistatic limit when sq = s/. From Eq. (|l^), it fol- 
lows that T c decreases linearly with a close to the QCP 
with a slope dT c /da T C (T /T C ) 2 / 3 for Ising spins. 

(ii) For Heisenberg spins, (3 = 2. T c o decreases with the 
reduction of quasiparticle coherence on approaching the 
QCP. T c is further reduced from T c q due to the increas- 
ing scattering rate. However, in contrast to the case of 
magnetic impurities where T c can be suppressed to zero 
at a finite concentration, we find that T c remains finite 
at the QCP as a consequence of the T-dependent scat- 
tering rate. From Eqs. (@,|2h at a = 0, it follows that 
f(T c )/T c = /3s;c"ln(A/T c )/s [l-|-c'ln(A/T c )] tends to a 
constant p = c"(3si/2sonc' of order unity for T c <C A, 
leading to T c ~ T cQ exp -[i/>(l/2 + p) - i/>(l/2)]. 
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FIG. 1. T c versus a for g 2 = 0.3 (lines) and 0.15 (tri- 
angles). T c o and T c obtained from Eq. ([17j) are shown for 
comparison. The peak in T c versus a scales with T*. 

The results obtained from numerical solution of 
Eq. (|l5|) are shown in Fig. 1. On the PM side, they 
are in agreement with those of Roussev and Millis ||. 
We next analyze how T c varies with a close to the QCP 
in the Heisenberg case. Eq. ( |l8| ) shows that T c q increases 
linearly with a, as the overall sign of the terms propor- 
tional to (T*/T c ) 2 / 3 has changed from the Ising case due 
to the real part of the self energy when (3 = 2. However, 
we find that this effect is subleading, and the dominant 
a-dependence of T c comes from the scattering rate in 
Eq. (|l7]) through the strong a-dependence of S" near the 
QCP, i.e. dT c /da oc -dS," '/da. From Eqs. (iP), we 
obtain, dT c /da ~ \/\[a, for y/a < (T c /A) ln(A/T c ); and 



dT c /da ~ 1/a, for (T c /A) ln(A/T c ) < y/E < (T c /A) 1 /3. 

In the FM phase, the spin rotation symmetry is broken. 
However, close to the transition, approximate Heisenberg 
symmetry is restored due to the small difference in the 
Fermi momenta which in turn, as discussed above, sup- 
presses the contributions of the long wavelength Gold- 
stone mode to the electron-spin fluctuation kernel. As a 
result, the superconducting phase boundary is approxi- 
mately symmetric near the magnetic QCP. Away from 
the QCP, the deviation from the Heisenberg symmetry 
leads effectively to a (3 value that is shifted downward 
from the Heisenberg value and a somewhat higher T c (see 
Fig. 1) on the FM side. Well inside the FM phase, the 
relative suppression of the fluctuation in the transverse 
channel makes the situation closer to the Ising case stud- 
ied above, resulting in a much higher T c 

To conclude, our results suggest whether a significant 
suppression of T c occurs near the QPC can be used to 
help identify the spin symmetry of the superconducting 
order parameter. We have shown that such a reduction 
occurs in the triplet case but is absent for singlet pairing, 
e.g. the s-wave pairing proposed in the weak ferromag- 
netic local Fermi liquid theory Q . In the singlet case, the 
spin fluctuations contributing to the self energy and the 
gap equation are identical and the dominant quasistatic 
pair-breaking effects in the quantum critical regime can- 
cel out as in the Ising case discussed above. While more 
experiments are clearly needed to determine the phase 
diagram near the FM QCP existing data in the AF 
case |]J show that T c indeed peaks near the magnetic 
QCP where pairing due to AF spin fluctuations is ex- 
pected to be of spin-singlet in nature. 
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